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Abstract 

Cq ; A shape ehange of the self-consistent mean-field induced by a configuration change is discussed 



X 



within the conventional constrained Hartree-Fock (CHF) theory. It is stressed that a single-particle 
level crossing dynamics should be treated carefully, because the shape of the mean-field in such a 



> 

finite many-body system as the nucleus strongly changes depending on its configuration. This situ- 
ation is clearly shown by applying an adiabatic assumption, where the most energetically favorable 
single-particle states are assumed to be occupied. The excited HF states and the continuously- 
connected potential energy curves are given by applying the configuration dictated CHF method. 
The effect of pairing correlation is discussed in the level crossing region. Triaxial deformed results 
in our Hartree-Fock-Bogoliubov (HFB) calculation with Gogny force nicely reproduce the available 



experimental data of Ge isotopes. From our numerical calculation, it is concluded that the CHFB 
state is more fragile than the CHF state in the level crossing region. 
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I. INTRODUCTION 



The HFB theory has been regarded to give a good starting point for a quantum description 
of an isolated many-fermion system like a nucleus. It defines an approximate ground state 
as well as a corresponding set of single-particle states, which have been used for the shell 
model calculations or the RPA treatments. Paying a price for breaking various conservation 
laws satisfied by the Hamiltonian, one obtains such an appropriate mean-field that might 
incorporate various correlations of the system as much as possible. Its many important 
theoretical concepts like a stability of the mean-field, a restoration dynamics of the broken 
symmetries, etc constitute an outstanding virtue of the self-consistent mean-field. 

Owing to intensive studies on the high-spin physics such as the band termination, band 
crossing, shape coexistence, identical bands and superdeformed bands, it has been shown 
that an introduction of many mean-fields to one nucleus successfully explains various excited 
states appearing near the yrast region. Namely, many mean-fields specified by an approx- 
imate angular momentum, aligned angular momentum, single-particle configuration and 
other approximate quantum numbers, provide us with physical, simple and satisfactory un- 
derstanding for a number of excited states better than a full quantum mechanical treatments. 



A rapid expansion of applications of the mean-field methods both with 



and without 



ion oi app 
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self-consistency gives rise to various important theoretical issues 
in the microscopic theory of many-body system. Among others, a property change in the 
single-particle states and its dynamical relation to the mean-field, a stability of the excited 
mean-field, a restoration dynamics of broken symmetries for the excited mean-fields and a 
dynamics of the single-particle level crossing are interesting subjects to be explored. 

In the level crossing region, there have been made many discussions on an applicabil- 
ity of the cranked mean-field theory Q], because there appears not only a large angular 
momentum fluctuation, but also a spurious interaction between two crossing orbits not at 
a given angular momentum but at a given angular frequency. In order to approximately 
restore the angular momentum conservation law by hand near the level crossing region, at 
an expense of the virtue of the mean-field theory stated above, it might be meaningful to 
eliminate the spurious interaction and to introduce a diabatic orbit 13[- With regards to 
the non-adiabatic (diabatic) effects due to the fast nuclear rotation, rotation-induced time- 

nn 

odd components of one-body density have been introduced and analysed in Refs. LL-HJ , 
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and detailed analysis [l|, |2|, |3j |4( shows their importance in the identical and superdeformed 
bands. A dynamical and quantum origin of these diabatic or non-adiabatic effects appeared 
in the mean-field theory also exhibits one of the most interesting problems in the theory of 
nuclear collective motion. 

Since these interesting problems have been mainly discussed qualitatively, it is necessary 
to study more deeply what actually happens in the self-consistent mean-field at the level 
crossing point, and how strongly it undergoes a change when it acquires additional defor- 
mation, additional angular momentum and when its single-particle configuration changes. 
Aiming to understand microscopically how the mean-field changes by itself as a smooth 
function of global variables, a new numerical method called the configuration dictated CHF 
method has been proposed in Ref. and recently applied for the rotational bands 16. 17^. 
In contrast with the conventional numerical method to solve the CHF equation, which are 
mainly applied on a selected set of mesh points on the global variables, our method intends 
to keep an identity of the CHF states and the calculation is made in a small step so that 
the resultant solutions are regarded as a continuous function of the global quantities like 
quadrupole deformation, angular frequency, etc. Applying this method, one may construct 
various configuration-dependent mean-fields as a function of rotational frequency or nuclear 
deformation, and explore what happens at the level crossing point. 

In this paper, we will pay special attention how strongly the self-consistent mean-field is 
affected by the configuration change, which is expected to occur at the level crossing region. 
This situation is well studied by comparing two different numerical calculations: one is the 
adiabatic mean-field theory and the other is the configuration dictated approach. By using 
the CHF and CHFB calculations the low-lying properties of Ge isotopes will be discussed in 
terms of the relation between high-?' intruder occupation and deformation, which has been 
Stud ,ed in the defonned m ea„-neld theory flHI^aid she., m „de. captions QQ 



in the A~80 mass region. 

This paper is organised as follows. Section |H] gives a brief outline of the theoretical 
method. In Sec. 1111} a specific feature of the adiabatic mean-field theory is discussed and 
low-lying excited states and continuously-connected potential energy curves are given by the 
configuration dictated method for Ge isotopes. The effect of pairing correlation near the 
level crossing region and a fragility of the CHF and CHFB mean-fields are discussed based 
on the numerical calculation. Triaxial deformed results in our HFB calculation with Gogny 
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force are compared with the available experimental data. Section HVl is devoted to discussion 
and summary. 



II. THEORETICAL METHOD 



In this section, we discuss the necessary formulae and methods for analysing structure 
changes near level-crossing region. For the sake of simplicity, we will discuss them within 
the constrained HF rather than the constrained HFB formalism. Inclusion of the pairing 
correlation is straightforward. 

In the present work, the self-consistent CHF(B) equation is solved using three dimensional 



23, 



24, 



25 



26|, El. The 



harmonic oscillator basis with the Gogny D1S interaction [2z 
Coulomb interaction and the center of mass correction up to the exchange terms are also 
taken into account. To save the CPU time in numerical calculation, we impose the Pe~ lnJz 
(z-simplex) and Pe~ m3y f (Sy) symmetries 28, 29], where P is the parity operator, e~ tnJt 
the rotation operator around i axis by angle n, and f the time reversal operator. Due to 
the z-simplex and ST symmetries, a mass asymmetry of the nucleus is allowed only along 
the x axis. To keep the center of mass motion fixed, we impose a quadratic constraint in 
the x-axis direction. The constrained HF equation to be solved is thus given as 



t{(H)+ l -C{(Q w )-tf+ l - 



a x (x)' 



0. 



with constraint 



(Q20) = Q, (x) = 0. 



(1) 



(2) 



A quadratic constraint of quadrupole moment is used to treat concave areas of the energy 
surface, /i is an input parameter which allows us to vary the expectation value (Q2o)- 
Meaning of C was discussed in Ref. jsO | - and is chosen to be 1.0 * 10~ 3 [MeV/fm 4 ] in our 
calculation. The Lagrange multiplier A is given by 



X = C(fi- (Q20)) = dE/dq, E=(H), 



(3) 



and an effective value of A is allowed to change during the iteration. In our numerical 
calculation, we take a x = 1.0 * 10~ 4 [MeV/fm 2 ]. 

For convenience, the constrained Hamiltonian is denoted by C = H — AQ20, hereafter. 
Equation is then expressed as (for simplicity, the center of mass correction constraint 



will be omitted in the following formula) 

5^(q)\Cmq)) = 0. (4) 

Hereafter an explicit deformation dependence of the CHF state \i[>(q)) is used instead of |). 
For \ip(q)), the constrained Hamiltonian is expressed as 

C(q) = H-X(q)Q 20 

= m)\C{qmq)) + Y.'M)t ] M)cM) -J2 e ^(l)cl(q)+ C{q) :, (5) 

fjj i 

where c 1 (q) and cj(g) are particle- and hole-creation operators satisfying 

5 M (gM?)> = 4(g)k%)> = 0. (6) 

Here and hereafter, the particle states are denoted by fi, v and hole states by i, j. If we 
do not distinguish, we use the letters k, I. The operator : C(q) : denotes the two-body 
residual interaction consisting of the normal-ordered product of four fermions in the two- 
body interaction with respect to \tp(q)). 

To understand the microscopic dynamics responsible for the origin of the quadrupole 
deformation, which will be discussed in the next section, the quadrupole operator is expressed 
as 

Q20 = q + ^{Q^cl&i + h.c.} + Qvv(si)c\c v - Y QiM)cjc\. (7) 

/ii /iv ij 

Here h.c. denotes the Hermitian-conjugation of the former term. The deformation- 
dependent particle-hole, hole-hole and particle-particle components Q^i(q), Qa{q) and 
Q^q) indicate how a microscopic structure of the self-consistent mean-field \ip(q)) changes, 
and what kind of single-particle configuration is favored as a function of nuclear deformation. 
The quadrupole deformation of the system is given by 

N 

g = £}<?«(?), (8) 

i=l 

N being the number of nucleons. 

In the present work, we solve CHF(B) equation by two ways. One is the conventional 
adiabatic method as will be discussed in Sec. IIII Al In this case, each iterative process is 
achieved by introducing a new density matrix constructed by the first N-lowest eigenvectors 
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(adiabatic assumption). A convergence is considered to be completed when the resultant 
density matrix is equivalent to the preceding one within a given accuracy. In this way, 
the most energetically favorable CHF(B) state satisfying a given constrained condition is 
obtained. In this method, one does not pay any attention to mutual relations between two 
neighbouring Slater-determinants with slightly different quadrupole deformation \if)(q)) and 
mq + Aq)}. 

To understand how the CHF(B) state undergoes a structure change depending on the 
quadrupole deformation, it is desirable to obtain \ip(q + Aq)) in such a way that it can be 
regarded as a smooth function of Aq. For this aim, we apply the second method called the 
configuration dictated method, which is briefly recapitulated below. Let \ip{q)) be a known 
CHF state satisfying a condition ('ip {q)\Q 20^ {q)) = <?• To find a new CHF state \ip(q + Aq)) 
which is supposed to be continuously connected with \ip(q)), we exploit a condition 

lim (<Pi(q)\(Pj(q + Aq)} = 8 id , i, j = 1, ■ ■ ■ , N, (9) 

Aq— >0 

where {\<Pi{q)),i = 1, • • • ,N} denotes a set of occupied wave functions constructing the 
single Slater determinant \ip(q)). That is to say, each small increment Aq is numerically 
adjusted by the maximum overlap criterion (JHJ) under a given accuracy, so as to maintain 
the identity of CHF states. In our calculations, Aq is so determined as to fulfil the condition 

|(^(g)|^(g + Ag))| 2 >0.9, i = l,.--,N. (10) 

In this way, the configuration specifying \if)(q)) is considered to be kept continuously as a 
function of q. Since the CHF state at q + Aq is dictated by the configuration of the preceding 
CHF state at q rather than the adiabatic assumption, this method is called the configuration 
dictated CHF method. 

This method is also generalised to get an excited HF state. Suppose there exists a HF 
state denoted by 1 0o) - O ne may then introduce a certain single Slater determinant \4> ), 
which is a simple np — nh state with respect to |0o) an d is called a reference state (RS). 
Here, it should be noted that the RS does not necessarily satisfy the HF condition. To get 
the HF solution |0 n p_ n h) locating at the nearest distance from \4> RS ), one may apply the 
following procedure in each iteration. 

1. Start the HF iterative procedure with the RS as an initial single Slater determinant 
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2. Suppose the nth single Slater determinant I0 n p_ nn ) is composed of a set of single-hole 
states {\ipf),i — 1, . . . , N}, where n denotes the number of iteration. Here, it should 
be noticed that {|v?™), i = 1, • • • , N} does not necessarily consist of the first N lowest 
states, but the occupied orbits. Using |</>™ n j 1 ), one may define the (n + l)th one 
body density P^p 1 ^; an d the (n + l)th single-particle Hamiltonian MP n j^)- 

3. With the aid of ^(P^p 1 ^); one ma y S e ^ a se t of single particle states {\ip™ +1 ),a = 
1, . . . , N, . . .}. In constructing the (n + l)th single Slater determinant 10^^)' one 
has to select i = 1, . . . , iV} out of {|<^a +1 ), a = 1, . . . , N, . . .}. To maintain a 
property of \(j) RS ) at each iteration, one has to find a new set of occupied orbits |<P™ +1 ) 
in such a way that there holds the following conditions 

l^f K +1 >i 2 >|(^f K +1 >| 2 , for aU k >j with j = l,2,..-,N. (11) 

4. Iteration of the second and third steps should be repeated till {|<p"), i = I, . . . , N, . . .} 
and {|<^™ +1 ), i = I, . . . , N, . . .} become equivalent within the required accuracy. When 
the iteration converges, one gets the excited HF state |0 n p_ nn ) whose configuration is 
dictated by that of \(j) RS ). 

The above-stated method is called a reference state (RS) method, and enables us to 
obtain many excited HF states. Applying the configuration dictated method by starting 
with various excited HF states, one may get many CHF lines (i.e. many potential energy 
surfaces as a smooth function of the deformation q) formed by continuously-connected CHF 
solutions, which will be discussed in Sec. 1111 Bl 

In the following calculation, both the adiabatic and configuration dictated CHF(B) are 
done in such a small step as the solution is considered to be a continuous function of the 
quadrupole deformation. Such a point-by-point heavy calculation is needed for discussing 
the dynamical change of nuclear system and the structure change of configuration-dependent 
mean-field as a continuous function of deformation. 

The single-particle wave functions are expanded in a three-dimensional harmonic oscil- 
lator basis up to the principal quantum number N Q = 8. The range parameters of Hermite 
polynomials have been optimised for each nucleus to reproduce the largest ground state bind- 
ing energy. The optimised range parameters thus obtained include some effects of higher 
major shells. Note that, no optimisation for each nucleus has been done in the most HF 



and HFB calculations [l|, |2j, La, |4( , although a larger configuration space is adopted. It is well 
known that the optimal parameters change depending on iVo, and then become rather stable 
as the shell number N increases j^jj]. In order to examine a size effect of the configuration 
space, low- lying spectra obtained by the RS method with N Q = 12 are compared to those 
with N Q = 8 for Ge isotopes, where the intruder orbit #9/2 plays an important role and pro- 
vides us with an excellent opportunity to study the shell structure and shape coexistence. 
The good agreement for excitation energies and deformations between Nq = 8 and 12 cases 
(see below) indicates a reliability of our discussion with N = 8. 

In our calculation, one set of range parameters has been used for each nucleus in numer- 
ically obtaining the CHF(B)-lines and excited states. A freedom of variation in the range 
parameters may lower the binding energy of CHF(B) states, which could lead to the effective 
mixing of different configurations. However, the use of fixed range parameters allows one 
to trace an evolution of the ground state configuration as a function of deformation, which 
makes the single-particle level crossing dynamics transparent. 

In our numerical calculation, a convergence condition is given as 

EK n) (?)-4 n ~%)|< 0.1 [i^n (12) 

k 

where is s.p. energy in the n-th iteration. 



III. NUMERICAL CALCULATIONS 
A. Adiabatic constrained HF 

The ground state for nucleus 70 Ge is obtained after an optimisation on triaxial deforma- 
tion parameters of the Hermite polynomials. Starting from the ground state, the adiabatic 
CHF calculations are carried out for nucleus 70 Ge shown in Fig. ^ As seen from Fig. L^a), 
the resultant quadrupole moment changes smoothly as a function of input quadrupole mo- 
ment parameter /z, except for three points denoted by A, B and C. The binding energy as a 
function of the quadrupole moment also exhibits drastic change at these three points, which 
is observed in Fig. L^b). 

To understand what happens in these three points from the microscopic point of view, 
the neutron and proton single particle energy levels are depicted in Fig. |21 as a function 
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of quadrupole moment q. Here Occ. and Unocc. are used to denote the occupied and 
unoccupied single particle orbits in adiabatic calculation. They are specified by the parity 
and signature quantum numbers (tt, a), because the asymptotic Nilsson quantum numbers 
are not good quantum numbers when the reflection symmetry is absent. For convenience, 
a subscript a, b and c for (ir, a) is used to identify the orbits responsible for the gaps A, 
B and C. In this mass region, gg/2 orbit plays an important role for the shell structure and 
nuclear deformation. In the following discussion, the configuration characterising the CHF 
state is labelled as \p1P2, ^1^2], where p t (%) denotes the number of protons (neutrons) in 
the oscillator shell Nq = 3, and P2 (77-2) the number of protons (neutrons) in gg/2 shell. 

The ground state configuration of 70 Ge is characterised by [12 0, 16 2], which is observed 
from Fig. |21 In Fig.^a), it is shown that the neutron unoccupied orbit (— , — ) a belonging to 
fa/2 or P1/2, and the occupied orbit (+, +) ffl belonging to gg/2 come close, as the quadrupole 
moment decreases from the ground state, and eventually cross with each other. Although 
p is equally discretized just before and after the gap A, the calculated quadrupole moment 
decreases dramatically. After the gap A, the orbit (— , — ) a becomes occupied and the orbit 
(+, +)o unoccupied in accordance with the adiabatic assumption. Namely, two neutrons 
are excited from g 9 / 2 to / 5 / 2 or pu 2 relative to the ground state, forming a new excited 
configuration [12 0, 18 0]. 

In other two gap regions, the same microscopic structure change takes place. At the point 
B, the neutron orbit (+,+){, belonging to gg/2 crosses with the orbit (—,—)& belonging to 
fa/2- After the crossing B, two neutrons occupying fa/ 2 move to gg/2, forming an energetically 
favorable configuration [12 0, 14 4]. Between B and C, the adiabatic CHF solution keeps 
its identity until the proton orbit (+, +) c belonging to g 9 / 2 , and (— , — ) c belonging to fa/ 2 or 
P3/2 cross with each other. After the crossing, another configuration [10 2, 14 4] is realized 
and the third gap C appears. 

A drastic change of the quadrupole moment just before and after the level crossing comes 
from three aspects. First is a smooth contribution on the quadrupole moment q due to a 
change of the parameter p, which is expressed by a straight line segment between A and B 
in Fig. [^a). Second is coming from a change of the configuration. Since the quadrupole 
moment of the system is produced by a sum of the quadrupole moment of the occupied orbits, 
and since the configuration change takes place between two orbits with different deformation 
character, i.e., one is usually of deformation driving and the other is of deformation anti- 
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driving, the configuration change between these orbits induces a substantial deformation 
change in such a finite many-body system as nucleus. Third is the contribution from the 
rearrangement of mean-field due to the configuration change. 

To demonstrate to what extent each contribution devotes to the gap, in the following, 
the above-stated three contributions are analysed quantitatively for the gap A. In this case, 
the change of quadrupole moment within an interval A/i=10.0 fm 2 is Ag=65.50 fm 2 , whose 
proton and neutron contributions are Ag p =20.08 fm 2 and Ag„=45.42 fm 2 , respectively. 
The gradient of the line segment between A and B shown in Fig. Q^a) i s estimated as 
follows; /j,b — /i J 4=360.0 fm 2 and qs — gA=184.7700 fm 2 . The proton contribution to the 
latter is (qs — <7a) p =81.3098 fm 2 and the neutron contribution is (g# — gA)n=103.4602 fm 2 . 
The gradient of the line segment between A and B is then given as: (qs — qAj/i^B _ 
/ia)=0. 51325, whose proton and neutron contributions are {qn — 9a)p/(a 4 b — A*a)=0. 22586 
and (qB — qA}n/{^B ~ A*a)=0- 28739, respectively. With the aid of the gradient between 
A and B, one may estimate the first contribution as Aq^ =2.2586 fm 2 and Aq^ =2.8739 
fm 2 within the interval Ajj, in gap A. In other words, the total quadrupole moment only 
decreases Aq^=Aqp + Agi 1 ' ) =5.1325 fm 2 when there happens no level crossing between 
occupied and unoccupied orbits. This smooth effect corresponds to 11.25% (~2.2586 fm 2 ) 
and 6.33% (~2.8739 fm 2 ) of Aq p and Aq n , respectively. 

To investigate the effect of the configuration change on quadrupole moment, the diagonal 
components Q mm (m can be hole or particle state) of the quadrupole operator are shown 
in Fig. 01 for the six specific orbits (four for neutron and two for proton) responsible for 
level crossings at the gaps A, B and C. A difference of the quadrupole moment between the 

occupied and unoccupied orbits at the gap A turned out to be Q(++) a — Q( ) a ~18.16 fm 2 , 

which is understood from Fig. Of a). Including the spin degeneracy, the second contribution 
is estimated as Aqi? ~ 36.32 fm 2 , which corresponds to 79.96% of the neutron effect Aq n at 
the gap A. On the other hand, Aqp =0.0 because there is no configuration change in the 
proton system. 

The remaining contribution, i.e., Ag„ =Ag n — Aq^ — Aq^ and Aqj? =Aq p — Aq^ — Aqf\ 
is regarded as an effect coming from the rearrangement of the mean-field, which should 
adjust by itself (self-consistently) in accordance with a change in the configuration. The 
third contribution turned out to be Agp 3 ' l ~17.82 fm 2 (88.75% of Aq p ) and Agi 3 ' ) ~6.23 fm 2 
(13.71% of Aq n ) for proton and neutron systems, respectively. Here, it should be mentioned 
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that the large rearrangement effect Aqp of the proton-field is induced by the neutron level 
crossing. 

Similar analysis are carried out for the gaps B and C. For the gap B where the cross- 
ing between neutron occupied orbit (— , — and unoccupied orbit (+, +)& takes place, the 
contributions from Aq^ and Aqf ) to Aq p are 10.30 fm 2 (37.00%) and 17.54 fm 2 (63.00%), 
respectively. The contributions to neutron field Aq n from Aq^ , Aq^f 1 and Aq^ are 13.1 
fm 2 (17.01%), 60.58 fm 2 (78.59%) and 3.4 fm 2 (4.40%), respectively. For gap C, the con- 
tributions to Aq v are 1.1 fm 2 (2.50%), 28.1 fm 2 (62.23%) and 15.9 fm 2 (35.27%), and the 
contribution to neutron field Aq n from Aq^ and Aq^ are 1.44 fm 2 (11.46%) and 11.12 
fm 2 (88.54%). From the above numerical analysis, it is clear that a large amount of shape 
change Ag^ 2 ) + Ag( 3 ) is induced by the configuration change. Namely, the shape of the 
mean-field in a finite system like the nucleus is strongly affected by a few orbits involved in 
the configuration change, and the self-consistency between the single particle wave functions 
and the nuclear mean-field should be taken into account appropriately. When one intends 
to discuss a relation between the diabatic and adiabatic orbits, one should take account of 
the above discussed large gap effects. The large gap appears not only in the nucleus 70 Ge, 
but generally occurs when there happens a configuration change due to the single-particle 
level crossing. 

B. Configuration dictated constrained HF 

In the previous subsection, it turns out that the adiabatic potential energy curve for 70 Ge 
consists of four different configurations, but has only one local minimum. By applying the 
configuration dictated method, we separate the adiabatic potential energy surface (PES) 
into several configuration- dependent curves with different low- lying minima. Fig. H] presents 
our results of the CHF potential energy curves for 70 Ge. Here filled circles are used for 
the adiabatic curves which just coincides with that in Fig. H^b), while crosses are for the 
configuration dictated curves. It is clearly shown that some excited minima are obtained by 
the configuration dictated approach, whereas only one ground state is in adiabatic PES. 

The procedure to obtain the configuration dictated curves in Fig. 0] are summarised as 
follows. First, the ground state of a given nucleus is obtained by optimising the three defor- 
mation parameters. Second, one constructs an appropriate np-nh state composed of some 
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occupied and unoccupied orbits. Taking the np-nh state as an initial trial wave function for 
the iterative process, one may obtain an excited HF state whose configuration is dictated by 
RS. Third, starting from the excited HF state, the CHF equation is solved under the max- 
imum overlap criterion by keeping the identity of CHF states. The continuously-connected 
potential energy curve is obtained for a given configuration. 

The ground state of 70 Ge has the configuration [12 0, 16 2], where the last two neutrons 
occupy one of g 9 / 2 orbits. After obtaining the ground state, a neutron 2p-2h state composed 
of an excitation from the most high-lying occupied orbit (+, +) a to the second low-lying 
unoccupied orbit (— , — ) a is constructed. Taking the neutron 2p-2h state as an initial trial 
wave function and applying the RS method, one may get an excited HF state with excitation 
energy ~1.1 MeV and the configuration [12 0, 18 0]. Starting from the excited HF state, 
the configuration-dependent curve is obtained satisfying the condition (jlUj) . a part of which 
just corresponds to the adiabatic solutions after the gap A shown in Fig. HJT>). 

Applying the same procedure, we obtain another excited HF state with excitation energy 
~ 2.7 MeV and the configuration [12 0, 14 4], where two neutrons are excited from the 
second high-lying occupied orbit (— , — )& to the most low-lying unoccupied orbits (+, +)&. 
The excited PES corresponds to the adiabatic solutions between gap B and C. The excited 
HF state after the adiabatic gap C, where two protons are excited from the most high-lying 
occupied orbit (— , — ) c to the third low-lying unoccupied orbit (+, +) c and two neutrons are 
from (— , —)b to (+, is obtained with excitation energy ~ 3.6 MeV and the configuration 
[10 2, 14 4]. 

To demonstrate systematically the discontinuity property of the adiabatic PES and the 
continuous property of the configuration- dependent mean-fields as a function of deformation, 
similar calculations are performed for 72 Ge, 74 Ge and 76 Ge as shown in Fig. For nucleus 
72 Ge, the gap in the adiabatic PES is formed by the crossing between one of the occupied 
(79/2 orbits and the unoccupied pi/% orbit. The ground HF state is specified by configuration 
[12 0, 16 4] and excited HF state by [12 0, 18 2], respectively. It is observed that there 
are two gaps and one saddle point in the adiabatic PES of 74 Ge. The gaps are also due to 
the crossing between the intruder orbit gg/2 and orbit in N = 3 oscillator shell, the most 
active orbits being 1/ and 2p. The saddle point is simply owing to the repulsive interaction 
between orbits with the same symmetries, e.g., both belonging to gg/2- Similar saddle point 
is observed for 76 Ge, and no gap appears in the adiabatic PES. 
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For Ge isotopes within excitation energy ~4 MeV, the low-lying spectra obtained by the 
i Nn = 8 and 12 are listed in Tabled together with the available experimental 
34j . In TableHl the energy differences between the excited states and ground 
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data 
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state given by AE = [B.E. (ground) — B. E. (excited)] are listed. 

The configurations specifying the excited HF states are also listed in Tab. HI Here, it 
should be mentioned that the configurations [piPp, nin 2 ] specifying the HF states are the 
same for Ng = 12 and iVo = 8. As noted in Ref. [35], a triaxial calculation is necessary to 
obtain the low-lying minima in the mass region A=60~80, since some of them are 7- soft 
nuclei. From our numerical calculation and experimental data, the ground and excited states 
of Ge isotopes are considered to have a triaxial deformation. Table shows that the binding 
energy for iVo = 12 is about 2 MeV lower than those in N = 8, while the excitation energy 
and deformations are quite similar for Nq — 8 and 12. From the above discussion, the size of 
configuration space is rather sensitive for the total binding energy, but not for the excitation 
energy, quadrupole nor triaxial deformations. 

In this mass region, the spherical or normal deformation are characterised by a config- 
uration where the Ng = 3 oscillator shell is occupied, whereas the larger deformation by a 
configuration with the occupied N = 4 (gg/2 orbits) oscillator shell ^Q,!^. For example, 
the first excited state in 70 Ge is characterised by a neutron 2p-2h state where two neutrons 
in the deformation-driving orbit g 9 / 2 are excited to the orbit in N Q = 3 oscillator shell. 
Consequently, the first excited state has smaller quadrupole deformation /3 2 (0.17) than the 
ground state (0.24). The second excited state is two neutrons excitation from Ng = 3 shell 
to the intruder orbit gg/2, resulting in the relatively larger (3% deformation (0.35). Since 
the third (fourth) excited state has similar configuration with the ground (second excited) 
state, the third (fourth) excited state has the similar /3 2 deformation with the ground (sec- 
ond excited) state. Upon the second excited state, another two protons are excited to the 
intruder orbit gg/2, forming the fifth excited state with superdeformed character (/3 2 = 0.42). 
Deformation character of low-lying minima for 72 Ge, 74 Ge and 76 Ge shows the similar situ- 
ation with those in 70 Ge (see Table HJ). It is clear from Tab. H] that our calculation indicates 
some low-lying excited states in this mass region. An experimental study on these shape 
coexisting phenomena would be of special interest. 
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C. Effect of pairing correlation 

We have also performed the adiabatic CHFB calculation with Gogny force to include the 
pairing correlation into the mean-field self-consistently. Figure EKa) shows the lowest HFB 
potential energy surface (PES) for 70 Ge, and the properties of the HFB ground state for Ge 
isotopes are listed in Tab. |TJ From this table, one may observe the binding energy increases 
0.59 to 1.34 MeV after having included the pairing correlation into the mean-field. In so 
far as the and 7 deformation are concerned, however, our numerical calculation does 
not show any substantial differences between the HF and HFB states, which is recognised 
from Tab. |TJ Since our triaxial deformed results both in the HF and HFB calculations well 
reproduce the available experimental data from nucleus 70 Ge to 76 Ge, the pairing correlation 
is not decisive for the ground state properties except for the binding energy. 

To study the effects of the pairing on the quadrupole deformation more deeply, in Tab. ITTl 
numerical values of (Q20) and (Q22) f° r Ge isotopes are listed for both the HF and HFB 
calculations. From Tabs. |T] and m one may state that the value (Q20) changes strongly 
depending on the neutron number, whereas that of 02 does not. In this sense, 02 and 7 
may not be regarded to be good quantities to specify the shape of the nucleus in this mass 
region. A comparison between the HF and HFB calculations in Tab. ITTl tells that the values 
of (Q20) and (Q22) generally become small (favor the spherical shape) when one includes the 
pairing correlation into the mean-field, except for a case in the value (Q20) of 74 Ge where 
the microscopic shell structure plays a role. 

In order to study what happens in the level crossing region after having included the 
pairing correlation, the neutron single-particle energies of the canonical basis for 70 Ge in the 
adiabatic CHFB calculation are shown in Fig. Efb). From this figure, it is observed that the 
configuration change induced by the crossing between two specific orbits still occurs in the 
CHFB calculation, even though the gaps, i.e., the missing regions in the adiabatic PES are 
not so distinct in comparison with the case of CHF in Fig. A smallness of the missing 
region may be understood by the following reason: two configurations having different shape 
are mixed up by the pairing correlation, and the wf-factor introduced by the BCS theory 
makes a concept of the configuration obscure. However, it turned out to be very difficult 
to get the excited HFB states when one applies the configuration dictated CHFB method, 
unlike the CHF case shown in Fig. 0] (This is a reason why only the ground state properties 
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in HFB calculation are presented in Table H}. In the case of CHF, the gap appeared in 
Fig. ^ does not indicate any difficulty of the mean- field theory, but is caused by the artificial 
adiabatic assumption put by hand. Applying the configuration dictated CHF method, one 
usually finds a continuation of the PES in the gap region, and obtains many CHF states with 
different configurations as discussed in the previous subsection. On the other hand, in the 
case of CHFB, the missing region remains no matter how much efforts we made in applying 
the configuration dictated method. From our numerical calculations, we may deduce the 
following conclusions: (1) the CHFB state can not be well characterised by the configuration 
defined in the canonical basis, and (2) the lowest PES in the CHFB is more fragile than 
that of the CHF in the level crossing region. 

Here, it should be mentioned that the configuration dictated CHF calculation also meets 
a difficulty of poor convergence or even non-convergence, when a pair of avoided crossing 
orbits comes close with each other. From our analytic and numerical studies it turned 
out that the CHF iterations diverge when the quantum fluctuations coming from the two- 
body residual interaction and quadrupole deformation become comparable with an energy 
difference between two avoided crossing orbits in the CHF one-body potential. The above 
stated difficulty is related to an applicability of the mean-field theory, when some part of two- 
body interaction could not be approximated by one-body potential successfully, but always 
remains there during the iterations. After including the pairing correlations, this situation 
might become more complicated, because there are competitions not only between the ph 
type two-body residual interaction and the HF potential, between the pp type two-body 
residual interaction and the pairing potential, but also their cross effects. These interesting 
problem will be discussed elsewhere 

IV. CONCLUDING REMARKS 

In the present work, the adiabatic and configuration dictated CHF(B) calculations are 
performed for Ge isotopes, and the following three points are discussed. 

First: by using the adiabatic CHF calculation, it is shown that the shape of the CHF 
mean-field in a finite system is strongly affected by only a few orbits involved in the configu- 
ration change. This situation appears not for the specific nucleus, but generally occurs when 
there happens a configuration change in association with the single particle level crossing. 
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One should thus carefully take account of the differences between the mean-fields accompa- 
nied with the diabatic orbits and with the adiabatic orbits. 

Second: the configuration dictated CHF method for Ge isotopes gives the low-lying 
spectra as well as the low-lying PESs, which nicely reproduce the triaxial deformed character 
of Ge isotopes. As far as the ground state properties of Ge isotopes are concerned, the pairing 
correlations are not decisive except for the binding energies. Since our numerical calculation 
indicates some additional low-lying excited + states, an experimental study on these shape 
coexisting states in this mass region would be interesting. 

Last: based on our numerical calculations, the PES in the CHFB is expected to be more 
fragile than that in the CHF in the level crossing region. Further study on the effects of 
pairing correlation near the level crossing region is needed in terms of the competitions 
between the two-body residual interaction and the one-body HFB potential, and between 
the average quantities like the -u-u-factors and their quantum fluctuations. 
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FIG. 1: Adiabatic CHF for 70 Ge: (a) the calculated quadrupole moment as a function of input 
quadrupole moment parameter /x; (b) binding energy as a function of the calculated quadrupole 
moment. 
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FIG. 2: Neutron (a) and proton (b) single-particle energies near the Fermi surface for 70 Ge. Occ. 
and Unocc. stand for occupied and unoccupied orbits in adiabatic CHF calculation, respectively. 
(ir,a) denotes the parity and signature, and its subscripts a, b and c represent the orbits responsible 
for the gaps A, B and C. 
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FIG. 3: Neutron (a) and proton (b) diagonal matrix elements of single particle quadrupole moment 
for the specific crossing orbits responsible for the gaps A, B and C. The symbol is same as that in 

Fig.m 




-605 



100 



200 



300 400 
q [fm 2 ] 



500 600 



FIG. 4: Results of the CHF potential energy curves for 70 Ge. Filled circles are used for adiabatic 
configuration and crosses for configuration-dependent curves. 
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FIG. 5: Results of the CHF potential energy curves for (a) 72 Ge, (b) 74 Ge and (c) 76 Ge. Filled 
circles are used for adiabatic configuration and crosses for configuration-dependent curves with 
independent configurations. 
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FIG. 6: Adiabatic CHFB for 70 Ge: (a) binding energy as a function of the quadrupole moment; 
(b) neutron single-particle orbits in canonical basis. Occ. and Unocc. stand for occupied and 
unoccupied orbits, respectively. The corresponding orbits are denoted by the same symbols as 
CHF calculation. 
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TABLE I: The triaxial deformed low- lying spectra by RS method for the binding energy (BE), 
quadrupole deformation parameter (fh), triaxial deformation (7), and the difference between the 
ground state and excited state energies (AE 1 ). Results with Nq = 8 and 12 are given for comparison. 
The configuration [piP2> ^l^] for Nq = 8 is same to the case in Nq = 12 and hence are not listed 
here. The ground state properties are presented in HFB calculation. The experimental data for 
binding energy [311 ] . quadrupole and triaxial deformation parameters are also listed. 

The energy is in MeV and triaxial deformation is in degree. 
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TABLE II: 


Expectation values of quadrupole operators Q20 


and Q22 for the ground state of Ge 


isotopes in 


MT^ anrl l-TEl^ ralmlatinrm 

111 Oj_L1V_1 111 1 J KjCXiIXj Llld UlUllO . 


(Omi flnH (Ooo) arp in 
2\j / Clll - K - L 22 / 1Ai 


fm 2 . 




Nucleus 


HF 




HFB 






{Q20) 


{Q22) 


{Q20) 


{Q22) 


70 Ge 


298.636 


116.276 


278.566 


105.906 


72 Ge 


405.678 


103.241 


324.727 


101.721 


74 Ge 


328.358 


99.223 


336.873 


98.719 


76 Ge 


383.277 


103.348 


354.340 


93.343 



25 



